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Building integrated photovoltaic (BIPV) systems are a relevant application of photovoltaics. In countries 
belonging to the International Energy Agency countries, 24% of total installed PV power corresponds to BIPV 
systems. Electricity losses caused by shadows over the PV generator have a significant impact on the 
performance of BIPV systems, being the major source of electricity losses. This paper presents a methodology 
to estimate electricity produced by BIPV systems which incorporates a model for shading losses. The 
proposed methodology has been validated on a one year study with real data from two similar PV systems 
placed on the south façade of a building belonging to the Technical University of Madrid. This study has 
covered all weather conditions: clear, partially overcast and fully overcast sky. Results of this study are shown 
at different time scales, resulting that the errors committed by the best performing model are below 1% and 
3% in annual and daily electricity estimation. The use of models which account for the reduced performance 


at low irradiance levels also improves the estimation of generated electricity. 


© 2014 Elsevier Ltd. All rights reserved. 
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1. Introduction 


Over the past two decades the number of grid connected PV 
systems has experienced an important growth. According to the 
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International Energy Agency IEA |1], between 1992 and 2010 the 
average growing rate of installed PV power was 39% in its member 
countries. This growth has been due to the inherent benefits of 
renewable energies and supporting polices by public entities, 
mostly based on feed-in-tariffs. The total cumulative power 
installed in IEA countries at the end of 2010 was 35 GW. Of this 
capacity, nearly 24% was integrated in or associated to buildings. 
Despite progressive reductions of feed-in-tariffs in most coun- 
tries, it is expected for BIPV systems to keep their growth rate. 
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With rising prices of retail electricity and decreasing costs of PV, 
grid parity with commercial electricity will soon become a reality, 
especially in urban residential areas. This fact, together with new 
policies that enhance self-consumption (currently under discussion 
in Spain [2]) points that building integrated photovoltaic (BIPV) 
systems will become more relevant in the future. Also to be 
mentioned is the new European Energy Performance of Buildings 
Directive [3], where PV can definitely play a new role in the so- 
called “nearly zero-energy buildings” which are described in the 
directive as buildings with a “very high energy performance. [...] The 
nearly zero or very low amount of energy required should be covered to 
a very significant extent by energy from renewable sources, including 
energy from renewable sources produced on-site or nearby”. 

The increase in installed PV power demands for simple, 
accurate models able to predict the electricity expected from a 
system. These models should be confronted with actual data from 
existing installations to validate them. However, in the literature 
there is a lack of specific studies about performance of PV systems 
integrated in buildings [4-6]. As a consequence there is also a lack 
of models which describe system behaviour with reasonable 
accuracy. 

Location determines the performance of a PV system in two 
ways: climate and obstacles surrounding the system. Available 
irradiance and ambient temperature are set by climate while 
obstacles may reduce the effective irradiance over the PV gen- 
erator. Orientation and slope of the PV generator surface also affect 
the effective irradiance received. However, while energy loss 
due to small deviations on orientation and slope from optimal 
position have low sensitivities [7], the impact of shading on 
system performance is relevant. 

Shadows over BIPV systems are casted by nearby objects such 
as trees, utility poles, surrounding buildings or even the building 
itself on which the system is integrated. Regardless their nature, 
these obstacles are characterized by irregular shape and closeness 
to the PV generator. The irregular shape makes its modelling and 
incorporation into the obstacle profile a complex and tedious task. 
The closeness of obstacles to the PV generator means that an 
important fraction of the sky dome is occupied by these obstacles 
and shading has a high impact on BIPV systems. Therefore, 
a practical estimation method of BIPV systems performance 
should include an accurate methodology for evaluating shading 
losses. This paper presents a methodology to estimate the elec- 
tricity produced by a BIPV system based on the combination of a 
novel model for the estimation of shading losses with well-known 
validated models for estimation of electricity production. The 
performance of this methodology has been assessed with a case 
study: The hourly electricity produced over one year by two BIPV 
systems integrated in the same building, one heavily shaded while 
the other is relatively free of shadows, has been estimated using 
the methodology and compared to the actual, measured values. 

The structure of the remainder of the paper is as follows. Section 
2 describes the procedure to estimate system output, paying special 
attention to the model for shading losses. Section 3 describes a BIPV 
system used to validate this methodology. Results and comments on 
this validation are given in Section 4. Finally, conclusions on the 
results and methodology usage are given in Section 5. 


1.1. Background 


Any method for the estimation of shading losses is based on the 
combination of two steps that can be considered as two indepen- 
dent processes: 


— First, the irradiance over the surface of interest - the PV 
generator - is obtained. This step needs as input two sets of 
data: a series of annual irradiation, typically hourly values, and 


a description of the obstacles surrounding the PV generator. 
The comparison of obstacle profile and sunpath determines at 
which hours the generator is obscured by shadows. 

— On the second step the electricity that can be extracted from 
the PV generator is estimated. This is expressed in terms of a 
dimensionless factor called shading loss factor: Ls, which can 
be added to the various losses also existing in a PV system. 


On the past decades several methods developed to address 
shading losses have been published on the literature. Some of 
these methods include the two steps described above or are 
limited to only one of the steps. 

Alonso-Garcia [8,9], Di Pizza | 10], Karatepe | 11], Kawamura [12] 
and Rauschenbauch [13] provide theoretical models to simulate 
the behaviour of PV modules and generators in conditions of non- 
uniform radiation. However, they do not provide methods 
to obtain the distribution of shadows over the surface under 
study. The majority of these models, the only exception being 
Kawamura’s one, take into account the effect of bypass diodes. 
Also, three authors (Rauschenbauch, Kawamura and Alonso- 
Garcia) include the effect of the inverse I-V characteristic of the 
shaded cell to estimate the I-V overall curve of the module or 
generated affected by shadows. 

Other authors, on the contrary, provide models to estimate the 
effective irradiation, that is, after incorporating shading effects but 
these models do not estimate the electricity losses in which the PV 
system incurs due the effect fn inhomogeneous irradiation. Drif 
[14] proposes a model based on geometrical description of 
obstacles in which considers each component of solar radiation, 
beam and diffuse, separately. Drifs model has been validated 
experimentally on a PV system located at Jaen University (Spain). 
Niewenda [15] introduces a software tool, SOMBRERO that allows 
to obtain the irradiation over any surface in the presence of 
obstacles. This tool is not limited exclusively to PV systems and 
can be combined with other tools (SUNCODE, TRNSYS) to evaluate 
the impact of shadows on passive applications of solar energy. 

Finally, Woyte [16] and Oozeki [17] analyze shading losses from 
the perspective of electricity production of the overall PV system, 
focusing their analysis on electricity losses rather than irradiation 
losses. Woyte studies the impact of shadows on a PV generator 
according to its electrical configuration: generator formed by 
several strings connected in parallel to one central inverter, each 
string connected to an inverter or AC modules. For large obstacles 
whose shadow covers a significant section of the generator Woyte 
does not find differences in performance between the three 
configurations. Oozeki introduces a model to estimate electric 
losses caused by shadows, as well as other losses, based on the 
monitoring of the system during one year. Oozeki’s model great 
advantage is that it does not require a characterization of the 
obstacles. On the other hand, its main drawback is that it requires 
a time span of one year and it cannot be used to estimated shading 
losses during the design phase of a PV system. 


2. Estimation of shading losses 


Fig. 1 shows an outline of the procedure proposed to estimate 
PV system output. Shading losses are the most relevant factor on 
BIPV systems performance and, therefore, the model to estimate 
them is the core of the methodology shown in the figure. The 
model for the estimation of shading losses has been developed 
according to the following objectives: 


— Accuracy: obstacle profile, despite its complexity, must be 
characterized accurately by the methodology; this will result 
in better estimation of associated losses. 
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Methodology for shading losses estimation 
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Nomenclature: 
a: orientation 
p: slope 
Q(0): global irradiation over the horizontal 
Blo P): in-plane beam irradiation 
Dla p): in-plane diffuse irradiation 
Relo, P): in-plane ground reflected irradiation 
Ta: ambient temperature 


Methodology for output energy estimation 


i 


Eac 


Tc: cell temperature 
Lco: optical losses 
Lee: electric losses 
Lex: shading losses 
Ls: system losses 
Eac: AC energy 


Fig. 1. Flowchart of methodology for system output estimation. 


— Simplicity: the methodology must be flexible and simple 
enough to describe obstacles, however complex they are, using 
few and easy-to-understand information. 


The models in Fig. 1 are described in detail below. 


2.1. Solar radiation model 


The three components of solar irradiation over an arbitrarily 
oriented and tilted surface can be estimated if the global irradia- 
tion over the horizontal plane is known. For this purpose scheme 
proposed by Liu and Jordan [18] is used: beam and diffuse 
components of horizontal irradiation are estimated from global 
irradiation through the existing correlation between clearness 
index and diffuse fraction. The use of a local correlation, when 
available, is preferred over the use of a global correlation [7,19]. 

Beam and albedo irradiation components over an inclined 
surface, B{a,?) and R,(a,f) respectively, during a given time period 
t are estimated in a straightforward manner as shown in Eqs. 
(1) and (2), where B,(0) is beam component of horizontal irradia- 
tion, G,(0) is global component of horizontal irradiation, 0s is 
incidence angle, @,, is solar zenith angle, a is orientation angle 
(0° for south orientation and positive to the west), 2 is tilt angle 
(0° for horizontal surface) and p is surface albedo: 


B-(a, p) = B,(O)- max(O, cos @5)/ COS zs (1) 


Ri(a,B) = G0): (1— cos 2): p/2 (2) 


Perez model is used to estimate diffuse irradiation [19]. This 
model has been selected over similar models due to its better 
performance [7,20-22]. 

The model described so far needs global irradiation over a 
horizontal surface to estimate the components of irradiation over a 
tilted surface which is, in turn, necessary to estimate shading 
losses. However, in BIPV systems, if irradiation is measured it is 
always available on the same orientation and tilt as the PV 
generator. Therefore, a procedure to estimate horizontal irradia- 
tion from in-plane irradiation is necessary since there are no 
models in the literature for the correlation between diffuse and 


global components of irradiation over tilted surfaces. A novel 
approach to obtain values of horizontal irradiation based on 
selecting days from a typical meteorological year (TMY) from the 
PV system location is proposed. This approach has the advantage 
that it relies on data sets that characterize statistically the local 
climate conditions and are increasingly available [23,24]. Every 
day of measured tilted irradiation is replaced by a day from the 
TMY. It is not possible to replace an actual day with the same day - 
same date - from the TMY because the diffuse fraction of both days 
will not be the same. Instead, every day is replaced by the day 
from the TMY which is most similar to it. Similarity between days 
is estimated by the function in Eq. (3), fs. The first term inside the 
brackets, pruypara, is the correlation coefficient between both 
days; this term assures that the days from the TMY will have a 
similar energy distribution than the days from the recorded set. 
While correlation coefficient provides days with similar power 
profiles, it does not guarantee that these days have the same 
irradiation. For this purpose a comparison between daily energies 
is included as second term; Gg rmy(0) is daily irradiation of the day 
from the TMY, Gadata(a,f) is daily irradiation of the actual day and || 
|| represents absolute value function. The search inside the TMY is 
restricted to a window with a size of 2w+1 days, centred on the 
day to be replaced and extended w days ahead and behind of it. 
From within this interval the day that maximizes fs is selected: 


fs = (rmy para + [1 — lI(Ga rmy (0) — Ga data(@, B))/Gadatala, AI} /2 (3) 


2.2. Shading model 


The proposed methodology for shading losses estimation 
is based on the following concepts: 


— Irradiation components are considered separately. Global irra- 
diation is separated into its three components: beam, diffuse 
and ground-reflected. 

— Obstacles and PV generator are modelled as three dimensional 
objects. On BIPV systems the size and distance of the obstacle 
to the generator are frequently similar to the dimensions of the 
generator itself. Therefore, it cannot be assumed that the 
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Fig. 2. Grid defined over the PV generator surface for shading losses estimation and 
translation procedure of obstacle coordinates between points inside the PV 
generator. 


generator is constrained to a single point; instead, the gen- 
erator is considered as a 2 dimensional surface. In addition, the 
PV generator will not be evenly lighted when shaded: some 
parts will be covered by shadows while others will be illumi- 
nated directly by the Sun. In contrast, the much simpler 
assumption that the generator is dimensionless only allows 
for two possible states: fully illuminated or fully covered by 
shadows. 

— Shading losses are estimated using geometrical considerations. 
Although the effect of non-uniform lighting conditions on the 
electrical behaviour of a PV system is not linear with the 
surface fraction covered by shadows [8-13] the model devel- 
oped only takes into account the geometrical configuration of 
the PV generator, not its electrical configuration. The use of an 
electrical model to account for shading losses in a PV system 
requires, in addition to estimate parameters associated to the 
model, that physical and electrical configuration of both PV 
modules (number of cells, series—parallel configuration, size 
and number cells protected by a single bypass diode) and 
generator (series—parallel configuration, placement of modules 
along their long or short side). This situation can be compli- 
cated even further if PV modules are placed following an 
irregular pattern, as it is the case on some BIPV systems for 
constructive or aesthetic purposes. On the contrary, the geo- 
metric approach proposed provides a less complex shading 
model: the only information needed concerning the PV gen- 
erator is its physical dimensions, orientation angle (a) and tilt 
angle (£). 


As shown in Fig. 1, shading losses are estimated following 
a two-step process. Both steps are described below. 


2.2.1. First step: irradiation over PV generator surface 

First step goal is to estimate the incident irradiation on any 
point of the PV generator surface. This step starts by defining a 
grid of points over the PV generator. Obstacles and generator sizes 
are assumed to be of similar magnitude, meaning that every point 
of the generator “sees” a different obstacle profile. The points in 
the grid are selected so that they are coincident with the corners of 
the modules in the generator (as indicated by the grey dots in 
Fig. 2). For every point in the grid an associated obstacle profile is 
calculated. An obstacle profile is a two dimensional function that 
assigns to every point in the sky a binary value: 0 if that point is 
occupied by an obstacle, 1 otherwise. The range of this function 
comprises values of azimuth from — 180° to 180° and elevation 
from 0° to 90°. The obstacle profile is defined as discrete function 
so it can be easily computed by software. The resolution is set at 
0.5° for azimuth and 0.25° for elevation. The obstacle profile 
function is formally defined in Eq. (4), where w is azimuth and y 


is elevation. 


0, (y) > obstacle 


=[-— 180,180 = [0,90 
1 otherwise wl 11801. = 10,20) 


OP(y, y) = { 
(4) 


Fig. 2 shows the procedure to obtain the obstacle profile for 
every point in the grid. Initially, the coordinates of all relevant 
obstacle points relative to one point of the PV generator surface, G; 
in Fig. 2, are known. To obtain the obstacle profile seen from point 
Gz all points in the obstacle are translated from the coordinate 
system centred in G; to a new coordinate system whose origin is 
G2. The translation vector is the coordinates of G4 with respect to 
G2; vector G2G,. For example, the coordinates of point O4, which 
belongs to the obstacle, with respect to G2 (r2) are obtained by 
adding the translation vector to the coordinates of O, relative to 
Gy: 


fo=fPy +G 2G, (5) 


Once all points of the obstacle have been translated, their 
rectangular coordinates are converted into spherical coordinates: 
elevation (y) and azimuth (y). The repetition of this translation 
process for all the points in the grid provides a set of obstacle 
profiles. This set of obstacle profiles contains all the necessary 
information to know the distribution of shadows over the 
PV generator surface at any moment. 


2.2.2. Second step: estimation of shading loss factor 

The second step estimates losses due to shading on a PV 
generator from the set of obstacle profiles generated in the previous 
step. It starts by evaluating individually the effect of obstacles on 
every component of solar irradiation: beam B, albedo R, circumsolar 
diffuse Decs, isotropic diffuse Dj; and horizon band diffuse Dy. Each 
component has a different behaviour regarding obstacles, by evalu- 
ating them separately five partial shading loss factors are obtained. 


— Beam and circumsolar diffuse irradiation share in common 
that they come from the Sun position in the sky. Therefore, 
their respective shading loss factors, Isp and Lsycs, are 
estimated in the same way. The grid superimposed to the 
PV generator divides its surface into rectangular regions. The 
behaviour of each region is governed by the four points at its 
vertexes: if one region has two or more of its vertexes in the 
shaded part of the PV generator this region does not con- 
tribute to the electricity production of the generator. The 
partial shading loss factor for beam and circumsolar diffuse 
is then defined as 


{number of shaded regions} 
{total number of regions} 


Leip = LsH,cs = (6) 


The reason why a shaded region does not affect the behaviour 
of non-shaded regions is because it is assumed that bypass 
diodes effectively isolate the effect that shaded modules have 
on non-shaded modules of the PV generator. By identifying 
every module with a region it is assumed that only one bypass 
diode is installed per module that limit the electrical power 
losses to the ratio between the area protected by one single 
diode and total PV module area [25,26]. This is a conservative 
approach since modules usually have several bypass diodes 
installed. However, with this assumption it is not necessary to 
know how many bypass diodes a module has, thus reducing 
the complexity of shading effects calculation [27]. 

— Isotropic diffuse irradiation comes from all points of sky 
hemisphere. For every point in the grid its associated loss 
factor lsy;s is the ratio between the solid angle of the obstacles 
and the total solid angle of the sky hemisphere seen by the PV 
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generator. The subscript in lss denotes that this loss factor is 
defined for a single point not for the whole PV generator. It 
must be noted that the amount of sky seen by a PV generator 
depends on its tilt angle. It is considered for convenience that, 
in the estimation of the solid angle occupied by obstacles, the 
radial distance is equal to 1. Therefore, it is omitted in Eq. (7). 
The shading loss factor for the whole PV generator, Lsrs, is the 
average of all partial factors obtained from the PV generator 
points defined by the grid: 


(solid angle of obstacles) 
(solid angle of sky hemisphere) 


SySy{OP(w, y); sin y - Ay - Ay} 
2(n—f) 


lsr ıs 
(7) 


Lsr s = mean (lsr is) (8) 


— Horizon band diffuse and albedo, although originated from 
different physical phenomena, have the same behaviour: they 
both come from the horizon line. The partial loss factors (Ist; 
ls4r) for every point in the grid are defined as the fraction of 
the horizon seen by the PV generator. A PV generator always 
sees half of the horizon line, an arc of z radians centred on 
generator orientation, not depending on its inclination and 
orientation. As before, the loss factors, Lsr and Lsrr, are the 
averages of the corresponding partial loss factors: 


(angle of obstacles) Zp (OP, y) Ay}... 
(angle of half horizon) 1 pie =O 


(9) 


lsr H = lsr R 


[spy = mean (lsr 4); Lsr,R = Mean (lsr r) (10) 


Finally, the overall shading loss factor, Lsp, is obtained from the 
weighted combination of the loss factors estimated for every 
solar irradiation component: 


Lsr = (Be - Lsr g + Decs + Lse.co + Dus + Lsr 1s + Dex 
- Lsr y +Re + LsrrR)/G (11) 


2.3. Power models 


The electricity produced by a PV system can be determined by 
estimating first the maximum DC electricity that can be delivered 
by the PV generator and later subtracting losses related to different 
physical phenomena: 


— Shading losses Lsp: these losses are described in Section 3.2. 

— Optical losses Lco: losses caused by the reflection of incident 
light when incidence angle is different than perpendicular. The 
model proposed by Martin has been used [28]. 

— Electrical losses Læ: electrical losses include losses due to 
voltage drops in wires on the DC side and PV modules 
mismatch. In this study electrical losses are included in the 
models for maximum DC electricity estimation and they are not 
considered separately. 

— System losses Ls: these losses account for the instantaneous 
DC-AC conversion efficiency of the inverter. The widely used 
model proposed by Schmidt has been selected [29] because its 
theoretical—practical approach, which allows to estimate mod- 
el’s parameters from actual measurements of inverter opera- 
tion. This model is presented in the following equation: 


N= Do/Pi = Po/(Po + losses) = Po /(Po +ko +kiı :po+k2 -p (12) 


where 7 is inverter efficiency, po is normalized inverter output 
power, p; is normalized inverter input power and k; are 


coefficients associated to conversion losses: kọ represents 
inverter self consumption, kı represents inverter losses whose 
dependence on power is linear (diode consumption, switching 
devices) and kz represents losses whose dependence on power 
is quadratic (wires, resistors, coils, etc.). 

— Thermal losses: since the effect of temperature on system yield 
is included in the models for maximum DC electricity estima- 
tion thermal losses are not accounted independently. 


AC electricity produced by a PV system is estimated by the 
following equation, which incorporates the losses previously 
defined: 


Eac = Epc.max(G, Tc) - (1 — Lsp) - (1 — Lco) - (1 — Lee) - (1 — Ls) (13) 


Maximum DC electricity produced by a PV generator, Epc max, 
depends only on operation conditions: irradiance and cell tem- 
perature. In this study two different models for the estimation of 
maximum DC electricity have been tested. The first model was 
proposed by Osterwald [30]. This model combines simplicity and 
accuracy and performs better than more complex models for 
crystalline silicon modules [31]. Osterwald model is presented in 
the following equation: 


Epc max = Pro: Gt - [1 — yp + (Tc —Tc*)]/C* a4 


where G; is irradiation, G* is irradiance at standard test 
conditions (1000 W/m7), Tc is cell temperature, Te is cell tem- 
perature at standard test conditions (25 °C), Pac is generator peak 
power at standard test conditions and yp is temperature coefficient 
of maximum DC power. 

Osterwald model has been used in combination with the 
nominal operation cell temperature model NOCT [32], presented 
in the following equation: 


Tc = T4 +(NOCT— 20) - G;/800 (15) 


where T4 is ambient temperature and NOCT is nominal opera- 
tion cell temperature. The model based on NOCT parameter has 
been selected because its performance is similar to more complex 
models for estimating cell temperature [33]. In addition, this 
parameter is available from the PV modules datasheets. 

During this study, Osterwald model was found to be inaccurate 
for low levels of irradiance. In consequence, a model that defines 
PV module efficiency as a function of irradiance was included in 
the study. The expression of Eq. (16), first proposed by Beyer in 
[34], has been selected, where e; are fitting parameters that must 
be evaluated for every PV system. The effect of cell temperature is 
the same proposed by Osterwald. 


n= Gr, Tc) = [e1 +€2 - G+ 63 - log(G,)]- [1-7 - (Tc—Tc*)] (16) 


3. Experimental validation 
3.1. PV system description 


The methodology proposed in Section 2 has been validated 
using a BIPV system located at the Technical University of Madrid 
(“Lopez Araujo” building, ETSI Telecomunicacion). The modules 
are placed above the windows of the building south façade acting 
as sunshades; in winter, when the sun is low on the sky, the 
sunlight shines through the windows helping to warm up the 
building; in summer, with the sun higher in the sky, the PV 
generator blocks direct sunlight into the building reducing thermal 
load. Fig. 3 shows this system. 

The system is formed by 282 PV modules, oriented 8° East and 
tilted 32°, distributed in 13 strings of serial connected modules. 
Nine of these strings are connected in parallel to a central inverter 
and the four remaining strings are connected to string inverters. 
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Fig. 3. PV system located on “Edificio López Araujo” building at ETSI. Telecomu- 
nicacion Madrid. (For interpretation of the references to colour in this figure 
caption, the reader is referred to the web version of this paper.) 


The system also includes a meteorological station with the 
following sensors: two reference PV sensors of the same technol- 
ogy as the PV generator (crystalline silicon) for the measurement 
of irradiance [35]! and cell temperature [36] and a PT-100 resistor 
for the measurement of ambient temperature. Both PV sensors are 
of same technology as the modules, with same orientation and tilt 
angles as the modules and located on the roof top on a place free 
of shadows. 

For this study information from two of the individual strings, 
S1 and S2, composed of 26 series connected PV modules each 
(outlined in red in Fig. 3), as well as data from the meteorological 
station, has been collected and analyzed. The reason for selecting 
these strings is that they are representative of BIPV systems where 
the PV generators act as sun shading elements, an architectural 
function that often implies partial shading in the PV modules. In 
fact, PV modules of string S2 are shaded during afternoon hours by 
the upper row S1. Thus, obstacle profiles for strings S1 and S2 are 
quite different, as shown in Fig. 4. These obstacles profiles have 
been traced from the most affected points of both PV generators: 
middle point of lower border for S1 and middle point of upper 
border for S2. 

These two strings have been chosen because of their similarity: 
same number of PV modules and connection to the same inverter 
model. The PV generators’ peak power at standard test conditions 
was characterized at the beginning of this study according to the 
procedure described by Caamaño et al. [37]. Since the strings are 
connected to the same inverter model, the MPPT strategy in both 
cases is the same. It is important to note that the MPPT efficiency 
of the inverters (and its associated error) is still present in the 
results presented in the next section. The conversion efficiency of 
both inverters was also characterized experimentally by measur- 
ing simultaneously inverters’ input and output during a day with 
clear sky conditions, therefore, the characteristic parameters 


1 It must be noted that the international standard describing the requirements 
for the measurement of irradiance has changed since the deployment of this 
system and that the irradiance sensor does not follow the requirements of the new 
version of this standard. In particular, the voltage of the cell used as sensor, which is 
polarized by a shunt resistor that allows measuring the current, sometimes exceeds 
the limit of 3% of Voc as indicated by the new version of the standard. This leads to 
a slight underestimation at high values of incoming irradiance. However, this effect 
is small in crystalline silicon solar cells. This is corroborated by the fact that the 
errors of this methodology at high irradiance are lower than the errors at low 
irradiance (see Section 4). 


obtained reflect inverters’ real operation conditions (PV generator 
power voltage curve). For this experiment a precision wattmeter 
(model LMG 500 by ZIMMER Electronic Systems GmbH) was used. 
The most relevant parameters of PV systems S1 and S2 are 
summarized in Table 1. The values in table were obtained directly 
from the documentation provided by the manufacturers and the 
installer except those marked with a star (*) which were measured 
by the authors. 


3.2. Data set 


The monitoring system records every minute instantaneous 
values of the following parameters: in-plane irradiance, cell 
temperature, ambient temperature, DC power and AC power. 

As mentioned before, the reference cell used to measure 
irradiance has the same orientation and slope as the PV genera- 
tors. As the radiation model described in Section 3.1 calculates 
irradiation over a tilted surface from horizontal irradiation it is 
necessary to apply the criterion from Eq. (3) to extract values of 
horizontal irradiation from a typical meteorological year of Madrid 
[7]. Different values for the time interval w have been tested: 0, 1, 
5, 10, 15 and 30 days. 

The power on the AC and DC sides of the PV system has been 
measured by the built-in monitoring systems of the inverters. In 
order to determine its precision, measured power data were 
compared with measurements recorded by a precision wattmeter. 
For both inverters the mean absolute error and its associate 
standard deviation are given in Table 2. 

The PV system under study has been monitored during one 
year: from March 2008 to February 2009. Before analysing model 
performance all records were checked in order to detect faulty 
data. Five days were excluded from the analysis because of 
partially or completely missing data. Additionally, 11 days were 
also excluded from the analysis in PV system S1 because the 
inverter monitoring system was not operating properly. Total 
number of excluded days amounts to less than 5% of the complete 
data set 


4. Results and discussion 


The methodology described in Section 2 has been used to 
estimate the electricity produced by the PV systems described in 
Section 3. The instant values recorded by the monitoring system 
were processed in order to obtain average values every 10 min. 
From these averages energy produced in spaces of 10 min has been 
calculated. The choice of a 10min interval is a compromise 
between precision and interval size [38,39]. 

Ten minute estimates and measured values were in turn 
averaged to get hourly values. A further filtering process was 
applied in order to remove values corresponding to in-plane 
irradiation lower than 5 Wh/m?. This filtering avoids the high 
uncertainty related to low values of power but it was kept as low 
as possible so that validation of the proposed methodology 
includes all weather conditions. 

Comparison between the results of the proposed methodology 
and real values has been carried out at hourly, daily and yearly 
scales. The accuracy of daily and hourly electricity production 
estimates has been determined with the symmetric mean absolute 
percent error SMAPE [40]: 


SMAPE(%) = (1/n) - 100 - Xillxi = yill/(Xi+Yi) (17) 


where n is the total number of points in the data set, days or 
hours; x; is the measured value and y; is the estimated value. 

The SMAPE parameter overcomes the problem found at early 
morning or late afternoon, when irradiation values are low and a 
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Fig. 4. Obstacle profiles for string S1, in black, and string S2, in light grey. 
Horizontal axis is azimuth and vertical axis is elevation. 


Table 1 
Relevant parameters for systems S1 and S2. 


Parameter name System 1 System 2 
Generator 
Slope 34 34° 
Orientation —8 —8 
Peak power* 2376 W 2263 W 
Modules in series 26 26 
Modules in parallel 1 1 
Module 
NOCT 47°C 47°C 
B —2.3 mv/°C —2.3 mv/°C 
4 0.004 °C7' 0.004 °C7! 
Cells in series 33 33 
Cells in parallel 2 2 
Inverter 
Praxi 2500 W 2500 W 
ko* 1.6% 1.2% 
k1* 3.0% 2.8% 
k2* 1.4% 2.3% 

Table 2 


Mean absolute errors and standard deviation of inverter built-in monitoring 
system. 


MAPE (%) o (%) 
Inverter S1 15 3.4 
Inverter S2 23 35 


small error of few watts in electricity estimation is overweighted 
by a more traditional accuracy estimator like the mean absolute 
percentage error (MAPE). Another advantage of SMAPE is that it is 
less sensitive than MAPE to errors in measured values (x;). For 
annual estimates of electricity production, since they are single 
values, the percent error 6 has been used: 


5 = 100 - (y; —X)/%; (18) 


Additionally, an external tool has been included in the study 
as a reference for a better assessment of the methodology 
performance: the commercial PVSyst software has been used. 
PVSyst has been chosen because it is a state of the art software 
widely used to analyze PV systems performance. It allows a 
thorough analysis of near shadings and PV strings can be split so 
that the behaviour of every module when shaded is modelled 
individually. 

Tables 3 and 4 summarize the performance in the estimation of 
AC power by the models described in Section 2 and for all values of 
w used to apply Eq. (3) as well as results provided by the reference 


Table 3 
Performance in the estimation of AC electricity for system S1 for different models 
and different window sizes. 


Model window size 6(Eaca) SMAPEd od SMAPEh oh 
(Eac) (days) (%) (%) (%)  (%) (%) 
Osterwald 0 —6.1 8.4 9.7 23.6 33.4 
Osterwald 1 —4.1 6.4 8.1 21.4 32.5 
Osterwald 5 —3.9 6.5 8.5 20.8 32.7 
Osterwald 10 -43 7.0 8.6 21.1 33.1 
Osterwald 15 —4.3 6.9 7.9 20.8 32.6 
Osterwald 30 —4.4 7.3 9.9 21.3 32.9 
Beyer 0 —0.9 6.9 7.9 23.5 35.3 
Beyer 1 —0.8 5.1 6.6 20.8 33.5 
Beyer 5 —0.5 3.3 4.1 19.1 32.8 
Beyer 10 —0.6 3.0 3.4 18.9 32.7 
Beyer 15 —0.6 2.8 2.8 18.8 32.6 
Beyer 30 —0.6 34 3.7 18.8 32.3 
Reference - —11.3 11.7 2.7 54.7 61.7 
Table 4 


Performance in the estimation of AC electricity for system S2 for different models 
and different window sizes. 


Model window size 6(Eaca) SMAPEd od SMAPEh oh 
(Eac) (days) (%) (%) (%)  (%) (%) 
Osterwald 0 —6.2 15.4 18.7 48.4 52.7 
Osterwald 1 -3.1 11.2 14.7 44.0 51.6 
Osterwald 5 —2.5 7.6 9.0 40.8 50.9 
Osterwald 10 —3.3 6.7 7.8 40.6 512 
Osterwald 15 —3.3 6.3 6.7 39.9 Sill 
Osterwald 30 —3.6 6.6 9.3 39.6 512 
Beyer 0 1.0 14.5 17.8 50.1 54.5 
Beyer 1 1.6 10.5 15.7 474 55.8 
Beyer 5 1.7 71 11.1 444 54.5 
Beyer 10 1.5 6.4 10.0 44.5 54.8 
Beyer 15 1.6 5.9 8.7 443 55.2 
Beyer 30 1.5 5:7 7.2 43.7 54.8 
Reference - —14.1 151 5.9 69.4 67.8 
50r 
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Fig. 5. Error in daily electricity estimations versus clearness index, comparison 
between Osterwald model and Beyer model for system S1. 


model (PVSyst tool). Table 3 shows the results for system S1 and 
Table 4 for system S2. For daily and hourly values, along with the 
average error (SMAPE) its associated standard deviation, ø, is also 
given. Production estimates for system S1 are more accurate than 
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Table 5 
Average error (SMAPE) in daily electricity estimation for system S1. 


Model All (%) Sunny (%) Partially overcast (%) Fully overcast (%) 
Osterwald 6.5 2.8 5.4 19.0 
Beyer 3.2 1.9 2.6 8.8 
Reference 11.7 12.3 11.6 10.2 
Table 6 


Average error (MAPE) in daily electricity estimation for system S2. 


Model All (%) Sunny (%) Partially overcast (%) Fully overcast (%) 
Osterwald 9.0 5.1 75 17.6 
Beyer 7A 45 7.6 16.6 
Reference 15.1 15.8 15.2 BS 
Table 7 


Annual irradiation and daily mean irradiation for every day type. 


Irradiation All Sunny Partially overcast Fully overcast 
(kWh/m?) (kWh/m?) (kWh/m?) (kWh/m?) 

Total 1641.5 814.2 732.6 82.4 

Daily (mean) 4.5 6.1 4.5 13 


those for system S2 in all cases. This result is not unexpected 
because system S2 has a more severe obstacle profile than S1. 

Comparison of results for different window sizes shows that a 
value of 0 for w, performs significantly worse than any other value 
for both systems; the remaining w values exhibit a similar 
performance with absolute differences smaller than 2% within 
models. For annual estimations a value of 5 days for w provides the 
best results in both systems. For daily and hourly values the lowest 
error corresponds to w equal to 5, 10 or 15 days, depending on 
model and system. The differences in performance for these three 
values are very close, with an absolute difference smaller than 
0.5%. However, the lowest standard deviation error is achieved 
with w=15 for all cases. These results show that windows 
corresponding to values of O and 1 for w are small in size and 
provide poorly matching days. A bigger value for w (5, 10 or 15) 
performs better since, for each day, a wider collection of days exist 
from which to select the most adequate day. Further increasing the 
value of w does not increase model performance anymore; a value 
of w =30 does not outperform smaller values, it only improves 
performance slightly for system S2. 

Regarding comparison between models, the Beyer model 
provides the best results for system S1 while for system S2 Beyer 
model performs better for annual and daily values while Oster- 
wald model performs better for hourly values. However, the 
differences between these two models are very narrow for both 
systems. When comparing with the reference model, both models 
tested in the study perform better. Interestingly, the reference 
model performs comparatively worse in the estimation of annual 
energy, the annual error being an order of magnitude larger. In 
addition, this model underestimates electricity production as 
annual error is negative for both systems. This fact is corroborated 
by the mean biased error (MBE), not shown on the tables, which is 
— 11.1% for system S1 and — 13.6% for S2. 

Fig. 5 shows the errors in daily electricity estimation versus 
clearness index for system S1 and Beyer and Osterwald models. 
Daily errors are bigger for small values of clearness index for 
Osterwald model while for Beyer model daily error is uncorrelated 
with respect to clearness index. These results corroborate the 


appropriateness to include a model that incorporates the effect of 
irradiance level on PV module efficiency. 

Finally, the performance of the proposed methodology has 
been compared with regard to meteorological conditions. The 
365 days analyzed were divided into three categories: 


— Clear days: the sky is free from clouds, diffuse fraction is below 
10% and the irradiance profile for the whole day has a very 
characteristic bell-like shape. 

— Partially casts days: the presence of scattered clouds in the sky 
covers and uncovers the sun repeatedly and irradiance changes 
its value very fast. 

— Fully overcast days: the sky is totally covered by clouds; 
irradiance is irregular over the day but with low values, never 
exceeding 600 W/m7?. 


Of the 365 analyzed days analyzed 136 were classified as sunny, 
164 as partially overcast and 62 as fully overcast. Three days could 
not be classified because of missing data. 

Average daily electricity estimation errors (SMAPE), separated 
into day types are shown in Tables 5 and 6 for Osterwald and 
Beyer models, using w=5 days as half window size. Average daily 
errors obtained from the reference model are also shown. The 
previous trends observed for all days also apply: electricity 
estimation is more accurate for system S1 than for S2 and the 
Beyer model performs better than Osterwald model. For these two 
models, performance for sunny and partially overcast days is 
significantly better than for fully overcast days, being performance 
on sunny days the best overall. The reference model, on the other 
hand, exhibits the opposite behaviour: it behaves better on fully 
overcast days than on sunny and partially overcast days. For fully 
overcast days, the reference model performs better than any other 
model for PV system S2 and it is comparable (slightly worse) to 
Beyer model for PV system S1. It must be noted that the impact of 
fully cast days in annual electricity production is significantly less 
than for clear and partially overcast days because of the differ- 
ences in irradiation for every day type. In this sense, annual global 
irradiation and its breaking down into day typologies are pre- 
sented in Table 7. Mean values of daily irradiation have been also 
included to enable comparisons between different weather con- 
ditions independent of the number of occurrences of every day 
during the year. Average daily irradiation for partially overcast 
days is 74% of irradiation available on sunny days while the 
percentage for fully overcast days is as low as 21%. When looking 
at accumulated irradiation during the year, 50% of total annual 
irradiation corresponds to sunny days, 45% to partially overcast 
days and only 5% to fully overcast days while 17% of days, 62 out of 
365, were classified as fully overcast. 


5. Conclusions 


A new methodology for estimating the electricity produced by 
BIPV systems has been presented. This methodology incorporates 
a novel model for shading losses estimation developed specifically 
for BIPV systems. This model for shading losses estimation is based 
solely on geometrical considerations, an electrical model for the 
PV generator or information about modules internal structure and 
wiring are not needed. The methodology is validated with models 
to estimate the useful electricity a PV system can produce. Two 
alternative models to estimate electricity production have been 
tested. Noteworthy, the model that estimates electricity without a 
direct measurement of PV generator peak power, not usually 
performed on BIPV systems, provides results comparable, better 
in some cases, than the model requiring an explicit measurement 
of such characteristic. 
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The methodology has been evaluated during one year with real 
data from two PV systems representative of BIPV. The evaluation 
has been performed over a wide range of irradiance values, 
starting from values as low as 5 W/m’, covering all weather 
conditions. The error obtained in annual electricity estimation is 
below 1% for the best performing model. The SMAPE of daily 
electricity is below 3% for the same model. It has been also found 
that including the effect of reduced performance at low irradiance 
levels improves the estimation of the electricity yielded by a PV 
system. When compared to a reference model widely used in PV 
system performance analysis that can be considered representa- 
tive of the actual state of the art, the methodology presented 
performed better than in annual, daily and hourly estimation of 
electricity production. 

The proposed methodology is simple, easy-to-use and can 
provide fast, accurate results at low costs since it does not require 
specific equipment. These features make this technology suitable 
for application in the supervision of BIPV systems. Although the 
methodology has been designed specifically for building inte- 
grated PV systems it can be applied also to other types of PV 
systems different than building integrated. 
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